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ABSTRACT 


Three cases of analytic data and one case of real data 
were numerically integrated using a 5-level baroclinic 
primitive equations model of the general Mug w I shgalleloe 
Experiments were performed using initial winds Gerived from 
the linear balance equation and also winds derived analyti- 
cally. The feasibility of using the linear balance 
equation to initialize the wind fieid was examined. In all 
cases, the forecasts remained meteorological and reasonably 
well-behaved. Nevertheless, the forecasts derived from 
initial winds generated by the linear balance equation 
excited large, operationally-undesirable inertial-gravity 
waves, while the forecasts from analytically determined 
initial winds remained virtually free of such small scale 
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Lion INGRODUCTION 


In ene “tela of operational numericah weather prediction, 
the tendency, in recent years, has been toward the develop- 
ment of sophisticated global prediction models. This has 
poem made possible by the rapid growth of computing capacity 
and developments connected with general circulation research. 

The purpose of this study was to examine a baroclinic 
Primitive equation model using a ghobal, staggered, spherical, 
Sigma coordinate system which could be used on an operational 
Gdeyeetomday basis by the United States Navy. Ineorder to do 
this, an input procedure was developed to allow a real time 
initialization from FNWC analysis fields. A set of analytic 
fields using an analytic spherical harmonic stream function, 
first presented by Neamtan (1946), was also developed for 
use as a controlled set of initial conditions. By using 
@hmese amalytic cases, errors in real data collection, 
edetysiSwama initialization, which are inevitable in pract: - 
cal meteorology, were avoided. The use of an analytic case 
also allowed the control of temperature and moisture distri- 
bution, predominant wave number, phase speed and wave 
amplitude. 

The objective was to isolate and correct problems in 
the model by using the controlled analytic invut with the 
ultimate aim of improving and extending the Navy's overall 


capability to predict weather phenomena on a global scale. 
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II. BAROCLINIC PRIMITIVE EQUATION MODEL 


The goverming differential equations, written in vector 
form, are similar to sets used by Smagorinsky et al. (1965), 
Arakawa et al. (1969) and Kesel and Winninghoff (1972). 
mae invegprations were carried out on a global, spherical, 
staggered, sigma. coordinate system using the conservative- 
type difference equations based on the work of Arakawa (1966). 
The complete set of finite difference equations are given 
in Kesel and Winninghoff (1972) for the non-staggered 
square grid. The Arakawa type spatial differencing was used 
to eliminate the spurious energy growth which can occur 
with the more conventional finite difference approximations 


to the nonlinear advection terms. 


A. PRIMITIVE EQUATIONS 
The vector equation of horizontal motion in the sigma 


(co) coordinate system for this model is: 
5 > > > an ace > > > > 
se (mV) + (VenV)V + 3 55 (OV) + f£(KXnV) + 1VO + oTaVn = F + D.. 
The Thermodynamic energy equation is: 
* j aa . 
ci + Ve(nTV) + 9 mG ) ~ TO 97 4 oV°Vn + to) = nH + D 
ot o0 Cp ah 
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The mass continuity equation is: 


T a 90 | 
eee ie 8 = 


Qs 


The moisture continuity equation is: 


3 ze Sonya 
se (mq) + V-(maV) + 7 55 699) = TQ + Dy 


These equations are supplemented by the hydrostatic 


equation 


Q 
Q- 


Q 
Q 


the equation of state 
G@ =,R7/ B 
and the dimensionless vertical coordinate 
© aa a a 
In the above equations, u is the zonal wind component; 
v, the meridional wind component; 7, the terrain pressure; 
co the vertical-velocity measure; T, the temperature; f, 


the Coriolis parameter; $6, the geopotential; q, the specific 


Wenreruy ssi, Cee diabatic neating; F, the frictional stress; 
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Di? Gwe stateral diffusion of momentum; D coe labperal 


ie 
errtrusion or heat; Das Vic lecerale airfrfusion ol molsture ; 
Q, the moisture source/sink term and P, the pressure. 
Furthermore, the formulations for frictional stress (7), 
lateral diffusion (D), heating (H) and moisture (Q) are 


identical to those described in Kesel and Winninghoff (1972). 


B. GRID 

The spatial finite differencing was performed on a 
staggered, spherical grid. The geopotential (¢), tempera- 
ture (T), specific humidity (q) and the vertical-velocity 
measure (w = - 6) were carried at the poles. The longitudinal 
and latitudinal grid increments were both five degrees. 
This gives 2520 points over the globe, with (¢,T,q,w) 
and (u,v) carried at 1260 points each (See Fig. 1). When 
¢,T,q or w variables are needed at (u,v) points (e.g. T 
in the friction term and w in the vertical advection of 
momentum term), they are defined as the average of the 4 


values at the surrounding (6,T,q,w) points. 


Cs VERTICAL LAYERING 

The present model divides the atmosphere into five 
Vayers, esesketched in Figure 2. The basic variables of 
the model are carried at the center of each ies. Addi- 
tional variables and conditions are specified at the 
interface between layers, as well as the upper and lower 


boundaries. 
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PEGURE 1. 
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M (South Pole) 


LOCATION OF VARIABLES 


i stands for ¢, T, q and w variables 
W stands for u and v variables 
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FIGURE 2. 


Level 


Computed Variables (Notational 
(at each level) 


Subscript ) 


UPPER BOUNDARY 


LOWER BOUNDARY 


VERTICAL LAYERING 
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O.4 
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hase 


In the above figure, sigma (0) is the dimensionless 


Vemereal coordinate, u, is the zonal wind component, 

yume thee meridionalswindecomponent, q is» the specific 
humidity, ¢ is the geopotential, tm is the terrain 
pressure and w is the vertical-velocity measure (-0c). 
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Inverraetons*erem-performed ain the Phillip's» (1957) 
sigma (c) coordinate system, where pressure, P, is normalized 
with the underlying terrain pressure, 7. The dimensionless 


vertical coordinate, sigma (0), is defined as 
o = P/t 


It foliows that 


fe) 1 at the earth's surface 


fl 


and fo} O at the top of the atmosphere. 


The zonal wind component (u), the meridional wind compo- 
nent (v), the temperature (T) and geopotential (¢) are 
Serriec at sigma levels .9, ./{, .5, -3 and .1. The specific 
humidity (q) is carried at sigma levels .9, .7 and .5. The 
vertical velocity measure, w = on is calculated diagnosti- 
early trem che continuity equation for the layer intertaces. 
The vertical-velocity measure (w) is carried at sigma lev ls 


=e .e, wemond .2. 


Die TLMESDIFFERENCING 

The two finite difference techniques that were used to 
step forward in time are the centered (leapfrog) difference 
scheme and the Matsuno (Euler backward) difference scheme. 
The finite difference equation for the centered technique is: 


al t-1l 


F = fF + 2A4t =— 


Lf 





Since the form has three time levels, it has both a physical 
mea Compucvational mode (Haltiner, 1971). This computational 
mode, which changes phase at every time step, causes the soiu- 
tions at adjacent time steps to become decoupled. Further- 
mere, cne™centered scheme is not feasible for the first time 
Svep. Therefore, the Matsuno forward time differencing scheme 
is used for the initial time step of each 6-hour integration 
cycle. Teen step not only reduces solution separation 
but also selectively dampens high frequency waves (Haltiner, 


1971). The finite difference equation for the Matsuno scheme 


ices 
t aFe 
Fe = FY’ + At yc 
ae ore 
EF = Fo + At at 


A time step of ten minutes was used for all experiments. 
This large time step is computationally stable when the 
Arakawa technique (Langlois and Kwok, 1969) of averaging 
quantities involved in the longitudinal derivatives is use l. 
Without this averaging technique, the von Neumann linear 
computational stability criterion (Haltiner, 1971) would 
require a 2.5 minute time step because the longitudinal grid 
distance at 85° north and south is only 47 kilometers. No 
averaging was done equatorward of 60° north and south. 

It should be noted at this time that lateral diffusion, 
friction, convective adjustment, convective condensation and 
large scale condensation are computed at every time step 


while heating is computed every 6 time steps (Kesel and 


Wearminghnoff, 1972). 
18 





Ce Hee eae OND LT TONS 


Pevae real and anasyute data were Wsed for initialfzation. 
The data was produced external to the main program on a 
Northern Hemispheric FNWC 63 x 63 grid. The analytic data 
with known phase propagation properties was produced 
following the work of Neamtan (1946), Gates (1962) and 
Heburn (1972). The real data was read from FNWC tapes. 
These fields were analyzed by the FNWC objective schemes. 

The main program was used to interpolate the data so 
that values at 5° latitude-longitude intersections were 
available over the Northern Hemisphere. The data was then 
reflected into the Southern Hemisphere. This procedure 
allows for an excellent test of the program. The model 
begins with fields which are mirror images of each other 
and they should remain so. In all cases, the input data 
consisted of the temperature analyses for the Northern 
Pemaesphere av 2 constant pressure levels distributed from 
1000 to 50 mb, height analyses at 10 of these levels, 
moisture analyses at \ levels from the surface to 500 mb, 
sea level pressure and sea surface temperature. Monthly 
mean surface temperature fields were used to derive an 
albedo field (Dickson and Posey, 1967). The terrain height, 
im Cm@e real data case, was that used by the current FNWC 
me@el. The terrain height, in the analytic data cases, wads 


set to zero. 


We 





imsemger to lige lWate sand evaluate certain numerical errors, 
it was decided to use analytic initial conditions for the 
mcorecam function ¥, frem which the true solution of the 
Pemplete vorticity equation in two dimensions could be 
Ge@ained. With the assumption of nondivergent horizontal 
flow, harmonic wave soiutions of the complete vorticity 
equation have been obtained for the sphere by Neamtan (1946). 
Geese solutions yield tne velocity of prepagation of the 


waves. The solution for the ww field was found to be 


} 


W = A sin(mi-vt) Pe enn 6) - Br°sine+ CP (sin 8) 


ee 


mrere A, Band C are constants to be determined, r is the 
radius of the earth, m is the hemispheric wave number, v/m 
is the angular phase speed of the wave disturbance 
(radians / second), an denotes the Legendre polynomial and 
ae represents the associated Legendre function. 

If we choose C = 0 (Gates, 1962) and n = m+1 [Buringt n 


and Torrance (1936), Kreyszig (1962)] equation (1) reduces to 


y = A sin(ma-vt )(2N!/25N!) sin 6 Coca” — Been Samed 
(2) 


Mae constant Ais arbitrary and proportional to the wave 
amplitude. But it was shown by Haurwitz (1940) that the 
solution obtained for ~W implies the existence of a velocity 
distribution over the sphere such that the angular velocity 


of the westerly current is made up of the sum cf three terms. 
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Te first term, B, is constant over the sphere; the second 
term varies along a meridian but is constant round a circle 
pieeravitude; and the third term represents a harmonic wave 
which is propagated zonally with a constant angular velocity, 


v/m, which is given by the formula 


7 n(n+1) - 2 22 
vin = Bat) ~ ACHAT) al 


where 2 is the earth's rotation speed. 
A reasonable meteorological pattern was obtained by 


sae =. mas 6 and v/m as -13°, 0° and 


choosing A as 1000 m 
+12° longitude per day respectively. The w field thus 
Semained was used to determine the™“initial fTlelds. 

The geopotential fields were derived by solution of the 
linear balance equation 


yoy = RE Sy EA F] (1. 


over the entire Northern Hemisphere. Equation (4) was 
solved using the ovale vaaetien iterative technique with a 
relaxation tolerance of one meter (Haltiner, 1971). 

The temperatures generated for the analytic cases were 
constant over each pressure surface and were consistent with 
the National Advisory Committee for Aeronautics (NACA) 
standard atmosphere (Haltiner and Martin, 1957). The eleva- 


tion 10 ,769 meters represents the tropopause in this 
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Eemescpnere. These foliowing relations were used in calculating 


the analytic temperatures. 


Tp(°K) = 288 - 0.0065 Z Z < 10,769 meters 

Beek) = 218 Z 2 10,769 meters 
7 P 0.19023 < 

Zp (m) = 44,308 [1 - (sO73735) ] Z 510,769 meters 

Zo(m) = 10,769 + 6381.6 1n(234-22) A 2a). 66) aie Se 


faethe foregoing equations, the pressure is in millibars and 
tme height Z, Bs in meters. In addition, the speeific 
humidity (q) at the lowest three levels was set to zero and 
the forecast procedure excluded friction, convection adjust- 
meme, large scale convective condensation, large scale 
eendensation, heating and terrain influences. 

The initial wind fields in the analytic cases were 
obtained in two ways. Winds obtained analytically from tl 2 
original yw fields wee called "analytic winds." Winds 
obtained from the linear balance equation were called 
"derived winds." It should be noted that "derived winds" 


were used for the real data case. 


A. ANALYTIC WINDS 
The initial wind fields for the "analytic wind" cases 


were derived directly by differentiation of the initial 


ae 





SS ——>  e_ —— = 


> Giese 








y field (Neamtan, 1946). The formulas for the non-divergent 


Wind components are 





—— uy {a sin(mi-vt )(-N)E (cos at Giese eye (cos 9)? 7 
2 NN 
- Br°cos o} (5) 
v= 1 A ( éNt) Sane Ccos 0) om cos (mA-vt ) (6) 
a NN 


Mmmere a is the earth's radius, Ais an arbitrary constant 
that is proportional to the wave amplitude, 8 is the latitude, 
v is the angular velocity of the wave, A is the longitude, 


mis the wave number and N = mtl. 


Ee SDERIVED WINDS 
The "derived wind" field was obtained from a wv field 


using the finite difference expressions 


S 

It 

{ 
9 |r 
>> 
ole 


and 


iss AY 
acos@ AdA 


where a is the earth's radius, A@ is the latitudinal 
distance increment and AA is the longitudinal distance 


increment. The w and ¢ fields were obtained in two parts 


ae 
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depending on latitude. Poleward of 25° north and south, the 
initiai* geopotential field, say ds» was retained and the ¥ 
field was obtained by solution of the linear balance 


equation 


Voy + Wwrve/f = 0°¢ /f 


Equatorwarc or 25° north and south, the field was deter- 


mined using 
hb = 
Y 1 ae 


where f is a mean coriolis parameter. Then a new geopoten- 
tial field (¢) was obtained by solution of the linear 


balance equation 


Zz 


Vo = (f£V°W + V-VEF) 


At 25° north and south, and $ were a combination of half 
the ~ and » from the poleward case and half the wv and 9 fr-m 
the equatorward case. It was hoped that this would reduce 
mice omeeruuee of the high frequency oscillations in the 


meopwes (Winninghoff, 1971). 


# 

By initial geopotential field is meant either the FNWC 
Objective amalysis of @ or a simulated ¢ field-obtained by 
solution of Equation (4) using the analytic function form 
of & Ciemation (1)). 
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IV. WAVE ANALYSIS METHOD 

A Fourier series was determined at eacn five degrees of 
Patatude around the latitude circle. The employment of such 
a series, known as harmonic analysis, is extensively used 
in the study of observational data (Jeffreys end Jeffreys, 
1956). With this technique, the phase angles and amplitudes 
@. Cech wave number around a latitude circle can be calcu- 
ated. A Fourier series can be expressed as follows 


(Heburn 1972): 


F(x) = A + i (A, cos mx + Be sin mx) 
= C. + Kee cos "nx — 6m) 
m 
where 
B A 
— m Z m 
m sin(ém) cos(ém) 
and | 


The first three experiments involved an input stream function 
of wave number six. The values of primary interest, there- 
fore, were 5 6 and Ce which are the phase angle and amplitude 
of wave number six. Other phase angles and amplitudes were 


extracted and examined, especially in the real data case. 
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Four experiments will be presented in this section. The 
first three of these experiments were performed with analytic 
fields. These analytic fields were derived from a stream 
function with a wave number of 6 and arbitrary constant, 


= The vertical temperature 


Moet equal to 1000 msec 
structure of these fields was determined by the NACA stan- 
dard atmosphere. Experiments I and II consisted of two 
36-hour forecasts each. The first forecast used "analytic" 
Winds and the second forecast "derived" winds. Experiment 
TII consisted of one 36-hour forecast using "analytic" 
Winds. Experiment IV was performed with FNWC analyses and 
consisted of one 36-hour forecast using "derived" winds. 
Experiment It. The analytic geopotential field used in this 
experiment was derived from a stream function with a wave 
number of 6 and phase speed of -13° longitude per day. Tre 
Mercial surface pressure analysis is shown on chart HE. Tie 
forecast surface pressure fields using both “analytic” and 
"derived" winds can be found on charts F-H. 

Experiment II. The analytic geopotential field used in 
this experiment was derived from a stream function with a 
wave number of 6 and phase speed of 0° longitude per day. 
The initial surface pressure analysis is shown on chart I. 


The forecast surface pressure fields using both "analytic" 


and “derived" winds can be found on charts J-L. 
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Table I compares the wave amplitudes and mean heights 
of the 36-hour forecast pressure field for experiment II 
(analytic winds) with those of the initial pressure field 
at selected latitudes. Tne amplitudes and mean heights in 
this experiment decreased at latitudes equatorward of 30° 
north (south) and increased poleward of 20° north (south). 
These amplitude variations are believed to be caused by 
memrrnear effects. 

In addition, Figures 6- 8 are plots of terrain pressure 
wersus’Torécast™nour at selected" latitudes. These terrain 
pressures are plotted every hour out to 36 hours. They 
Show small amplitude internal gravity waves with periods 
ranging from 6-12 hours. These oscillations are believed 
vo be due™to the geostrophic adjustment mechanism and™are 
peculiar to the type of analytic field used. 
Peverimenv'lll. Thewanalytic™ifield used™in thi's"experimet 
was derived from a stream function with a wave number of 
and a phase speed of 12° longitude per day. The initial 
surface pressure analysis is shown on chart M. The fore- 
weet <upeace pressure fields using "analytic" winds only 
can be found on charts N-P. 

Experiment IV. This experiment was performed using FNWC 
objective analyses for 00Z February 8, 1973. The initial 
surface pressure analysis is shown on chart A. The surface 
pressure analysis for 122 February 8, 1973 is shown on 


chart B. The 12-hour forecast made by the current FNWC 


ani, 





model is shown on chart C. The 12-hour forecast made by 
the global model is shown on chart D. 

it sneuld be noved™ that Experiment IV included terrain, 
frictional stress, heating and moisture effects as discussed 
in Kesel and Winninghoff (1972). 

Table II summarizes the results of all experiments 
performed. The analytic phase speeds are compared to the 
actual phase speeds for all experiments using analytic data. 
The table also includes A which is proportional to the 
mp litude of the: analytic wave, Bwhich is a function of 
the phase speed and wave number, the method of balancing 
and the time required to make a 36-hour forecast. 

The forecast fields in all the experiments using analytic 
data showed a considerable tilt backward at high latitudes 
in the phase propagation of the wave. This was to be ex- 
peeved Sinee the Arakawa averaging techniquewtends to smc7th 
mae pradients at high latitudes. In addition, the Arakaw. 
Veema2que™e.ves an effective Ax which is comparable To CIF ut 
at low latitudes, thus as the wavelength decreased toward 
the eens the phase speed also decreased (Gates, 1959). 

The result of this differential movement, which was 
more pronounced in the "derived" wind cases, was the forma- 
tion of closed highs and lows at the higher latitudes which 
propagated equatorward. This distortion was aggravated by 
nonlimear effects introduced after the field ceased to be 


hermmenic in the longitudinal direction. 
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Furthermore, in all the experiments using the "derived" 
Wane, Nigh frequency imértial gravity waves were generated 
Gue to the initial imbalance between the mass and wind 
fields in the tronics. These gravity waves caused a rapid 
deterioration of the harmonic six-wave pattern and added 
appreciably to the probiem at the poles. The shear which 
developed at 25° north (south) was attributed to the method 
used to blend the winds between the tropics and mid-latitudes. 
Miemuarper error dasthe phasewspecds calculated for the 
"derived" wind cases were also attributed to this initialization. 

Since these experiments were performed using a multi- 
level primitive equation model which allows divergence, 
Equation (3) was satisfied only approximately. Also, 
Rossby (1939) has shown that the presence of divergence in 
a barotropic atmosphere will slow the rate of wave propaga- 
tion, especially for small values of wave number. It was 
not surprising, therefore, that the actual phase speed we: 
always less than expected. 

Finally, a thorough examination was made of Experiment 
Ly og determine Mice ete hel Cue G YyoOtecne program, Char 
is,°ai the inertial-gravity motions were being controlled 
Pemiistically so that the model produces consistent meteor- 
Ological appearing results. Figures 3 - 5 rive the surface 
pressure oscillations at each hour of a 36-hour forecast 
Por selected latitudes using the full global model. In 
addition, Table III compares forecast pressures of the 


global model to both the verifying analysis and forecasts 
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made by the FNWC model for selected lows and highs. [In 
interpreting these results one should keep in mind the 
erude initialization. Aiso, for those not accustomed to 
seeing values of terrain pressure at each hour in a primitive 
equatton model, these oscillations, although larger than 
normal at these points, are by no means unusual for such a 
meeriyesinittalized integration. There is no artificial 
smoothing done and recall that the run begins with a linear 
meerance wind. in addition, there are none of the various 
devices which operational experience has shown is necessary 
Me five an acceptable short term product, i.e., not calling 
mecenvecvive adjustment for the first several time steps. 
Mevervheless, these results are Stable and retain the 


dominant meteorological scale without difficulty. 
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TABLE I 


Latitudes Amplitude of Wave #6 (Ce) Mean Pressure (Cy) 
Initial 36-Hour Forecast Initial 36-Hour Forecast 
Pressure Pressure Pressure Pressure 
Field Pale Field Field 
Equator 54 yy | 1051 1040 
10° N G2 Lae 14035 1034 
20° N 20). 1 24.0 GAS Oi 
30° a 302 Biel 2 995 995 
40° N 2d PAS 5, SES eH ET 
50° N Om. 130 953 956 
60° N aa | ome 935 937 
70° N 1 1.5 922 gah 
80° N 5 eal 5 Sie 914 919 


Note: Data is from Experiment II (Analytic Winds) 


TABLE I. Wave Amplitude and Mean Heights for Selected 
| Batitudes and Forecast Time. 
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Air» ih 


Location of Initial Surface Verifying FINWC Global 
Pressure System Pressure Analysis Surface Pressure 12-hour 12-hour 
en initial 00Z 8 Feb 1973 Analysis Surface Surface 
Surface Analysis 12Z 8 Feb 1973 Pressure Pressure 
Analysis Analysis 


65°N 135°W 1042 (H) 1041 1040 1042 
foe), 15°E 971(L) it D 976 a5 
SSeN 25°W 991(L) 990 989 989 
65°N 70°W Soon) 392 992 Yige 
He-N WE5°r 966(L) 962 959 962 
5 Saal, 105° E 1049(H) 1046 1043 WO 32 
Sool B45 °F 974(L) 976 vite ce 
poor 553W 1031 (H) 1030 1029 1036 
50°N 100°W Os ( (H) HOS J, Si 1040 
50°N 160°W 1004(L) 1009 1004 1005 
B5-N 175°U 1030(H) 1030 1029 1S) 5 iL 
Soon 20° 1006(L) 1009 1007 1003 
So°N soc HOTZ CL ) 1006 LOWry 1010 
TABLE III. Comparison of 12-Hour Forecast Pressures of 


the FNWC Model and the Global Model Against 
the Verifying Pressure Analysis (Experiment IV). 
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CHART B. Verifying Surface Pressure Analysis 
(Experiment IV) 
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l2-Hour Forecast Using FNWC Model 
(Experiment IV) 
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CHART D. l2-Hour Forecast Using Global Model 
(Experiment IV) 
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CHART F. 1l2-Hour Surface Pressure Forecast Using 
Analytic and Derived Winds (Experiment I) 
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CHART H. 36-Hour Surface Pressure Forecast Using 
Analytic and Derived Winds (Experiment I) 


U7 








CHART I. Initial Surface Pressure Analysis 
(Experiment IT) 
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CHART J. 12-Hour Surface Pressure Forecast Using 
Analytic and Derived Winds (Experiment IT) 
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CHART K. 24-Hour Surface Pressure Forecast Using 
Analytic and Derived Winds (Experiment IT) 
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CHART L. 36-Hour Surface Pressure Forecast Using 
Analytic and Derived Winds (Experiment (II) 
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CHART M. Initial Surface Pressure Analysis 
(Experiment IIT) 
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CHART N. 6-Hour and 12-Hour Surface Pressure Forecast 
Using Analytic Winds (Experiment IIT) | 
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CHART 0. 18-Hour and 24-Hour Surface Pressure Forecast 
Using Analytic Winds (Experiment IIT) 
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CHART P. 30-Hour and 36-Hour Surface Pressure Forecast 
Using Analytic Winds (Experiment III) 
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VI. CONCLUSIONS 


ites Casco Or analyvlewdava and ome case of realmdata 
were numerically integrated using a 5-level baroclinic 
primitive equation model of the general circulation. 
Experiments were performed using winds derived from the 
linear balance equation and winds derived analytically. 
The feasibility of using the linear balance equation to 
initialize the wind field was examined. In all cases, the 
forecasts remained meteorological and reasonably well-behaved. 
The forecasts using winds analytically derived were virtually 
free of small scale inertial-gravity motions,’ while the 
forecasts weir. winds from the linear balance equation 
excited intertial-gravity waves which were large and unde- 
Sirable for operational forecasts. 

Therefore, an important question to be answered is wi: it 
method of balancing should be used to initialize the glot 2l 
Mocel™ so thav it-willl not suffer during the early part o 
the forecast run from excitation of excessive inertial 
gravity motions. A number of solutions have been proposed 
by Myakoda and Moyer (1968), Nitta and Hovermale (1969), 
and Winninghoff (1971). 

The method examined by Winninghoff used the equations 
in the model itself in an iterative sense either at a fixed 
time level or even in a four dimensional sense in which 


data is assimilated into a running model. This allows for 
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the natural adjustment mechanism itself to achieve the 
desired balance. A time scheme such as the Euler-backward 
may be used for selective damping of high frequency waves. 
The iterative method proposed by Winninghoff has the advan- 
tage of mathematical simplicity and complete consistency 
with the prediction model. Unfortunately, Winninghoff 
(1971) estimates that the iterative procedure must continue 
for an equivalent of an 18 to 24 hour forecast. This, of 
course, is not operationally feasible at the present time. 
Consequently, the technique which is now being tested by 
FNWC utilizes the best vertical mass structure possible obtained 
by the WVavest variational analysis techniques available 
and then solves the complete balance equation at each level 
iimeaderm wo get the initial wimd fields. This initialization 
along with other helpful operationally tested devices rch 
Bowomeeuvluang, turning friction on slowly and not calling 
en cConmmectime adjustment for theefirst» several time steps 
Should be a satisfactory interim solution to this most 


Gitereult problem. 
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